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ABSTRACT 

This  paper  presents  ?n  algorithm  for  solving  interval  linear  programing 
problems.  The  algorithm  is  a finite  iterative  method,  which  in  each  iteration 
solves  a full  row  rank  interval  linear  programing  problem,  with  only  one 
additional  constraint.  The  solution  and/or  problem  chosen  appears  to  be 
conputationally  more  efficient  than  that  in  the  Ben- Israel  and  Robers  algorithm. 
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# 


Introduction 


The  general  linear  programing  problem  (LP)  may  be  defined  as: 

CLP) : Max  cTx 

subject  to 
Ax  * b 
x > 0 

where  the  matrix  A and  the  vectors  b and  c are  given. 

Every  (LP)  problem  is  equivalent  to  an  interval  linear  programming 
problem  (IP)  of  the  form: 

(IP):  Max  cTx 

subject  to 

b < Ax  £ b+ 

where 

*•(?)  b'-W 

and  U > 0 is"sufficiently  large"  (e.g.  non -Archimedean  transcendental  for  the 
unbounded  case). 

Because  of  the  above  reason  and  since  the  constraint  system  in  many 
linear  programming  problems  arise  in  (IP)  form  (wholly  or  partly),  special 
algorithms  for  (IP)  problems  are  important,  and  may  be  more  efficient  for  those 
particular  problems  than  general  algorithms  for  (LP)  problems. 

0.  Preliminaries  and  notation 
0 - the  empty  set 

(x)  - the  set  consisting  of  the  single  element  x 
{x:  f(x)}-the  set  of  x satisfying  f(-) 
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iff  - if  and  only  if 

Rn  - the  n dimensional  real  vector  space 

For  any  xf  y c Rn: 

x >_  y - denotes  xi  > y.  1 £ i £ n 
x ± y - denotes  .2  x^y^  ■ 0 

For  any  subspace  L of  Rn: 

C*  « {y:  xly,  for  all  x c L}  the  orthogonal  complement  of  L 
x+L  - the  manifold  (x  ♦ 1,  1 c L) 

R"1™  - the  space  of  m x n real  matrices 
K™*1  -(xtr;  rank  x - r) 

I - the  n x n id  ntity  matrix 
n 

e.  - the  i**1  column  of  I„ 
i n 

e ■ ? e. 

i-1  1 

For  any  A e R^01: 

T 

A - the  transpose  of  A 

R(A)  » {y  e Rn  : y ■ Ax  for  some  x e Rn}  the  range  of  A 
N(A)  * (x  e Rn  : Ax  ■ 0)  the  null  space  of  A 

I 

A * a right  inverse  of  the  matrix  A 

An  interval  linear  programning  (IP)  is  defined  as 
(IP):  Max  cTx 

subject  to 

b < Ax  < b+ 
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where  c = (Cj),  b"  - (b.‘),  b+  » (b*)t  A - (a^)  (1  < i < m;  1 < j < n) 

are  given  and  b*  ^ b+. 

Any  x e Rn  satisfying  (2)  is  called  a feasible  solution  of  (IP).  If  (IP)  has 

T 

feasible  solutions  and  {Max  c x:  x feasible)  is  finite  then  (IP)  is  bounded. 


Lemna  1 [1] : 


Lenina  2 [ ! ] : 


A feasible  (IP)  is  bounded  iff  c e N(A) 

If  A is  of  full  row  rank  then  the  optimal  solution  x* 
of  (IP)  is  given  by: 


x*  = Z b:  A*  ♦ Z b+  A*  + I!  lO.bT  + (1  - 0.)bj  A*  + N(A) 
i l 1 " 1 1 id*  1 1 i c i o 11  1 1 1 


where 


and 


r = { i:  C Aj  < 0} 

I+  = {1:  CTA*  > 0} 

1°  - { i : CTaJ  = 0} 

° < Qi  < i,  i e 1° 

If  i 

A^  * 1th  column  of  A 


In  this  paper  an  algorithm  for  solving  (IP)  problems  in  the  general  case  where 
A is  not  of  full  row  rank  will  be  described.  We  shall  first  solve  a special 
class  of  (IP)  problems  where  the  coefficient  matrix  can  be  separated  into 
A 


T . 


(V)  where  F is  a full  row  rank  matrix,  and  h1  is  a single  vector. 
This  result  forms  the  basis  for  our  algorithm  to  solve  the  general 
(IP)  problan,  with  coefficient  matrix  , where  H is  any  matrix.  See 

also  Ben- Israel  and  Robers  [6],  [7]. 


1.  - The  Subproblem 

In  order  to  ^olve  the  general  (IP)  problem  we  will  describe  in  section  2 
an  iterative  algorithm  which  will  in  each  step  make  use  of  the  method  described 
below  to  solve  subproblems  of  the  form: 


4- 


0) 


Max  cTx 
subject  to 

d~  < Fx  < d* 


(2) 

™ g < hTx  < g4 

where  F is  of  full  row  rank. 

Let  us  substitute 

(*)  * - Fx 

C1) » (2) , (3)  in  order  to  obtain 

(5)  Max  cTFfz 

subject  to 

C6)  d”  s z < d4 

(?)  f“  < h^z  < g4 

where  F1  is  the  right  inverse  of  the  matrix  F. 
Let 


(8) 

Thus  we  have 

(9) 


pl-i'i  ifhTF*»0 

d\  - if  hTFj  < 0 


0 < ^ < d[  - d[ 


If  fj  refers  to  indices  i for  which  h?Fj  > 0,  f0}  to  those  for 


which 


•i . f 


!\  ■'.>  (’f  to  indict*s  for  which  hIpj  < n»  then  (5)  can  be  written  equivalently  as 

m J0  eVA  ♦ <9  ♦ Z - ip  . 


(10)  Max  l e^'d l . I cTF'd*  . I cTF*2f.  ♦ H-1)cTf£. 


SO 

Let  us  denote  by 


n(?)  " I h Fi^i 4 di>  ♦ I bTpJ(dJ  - - hV, 

♦»0 


n(0)  - I hTFjdT  ♦ I hTFjdt 


♦.o 


CH) 

then 

(12) 
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Using  (10U11),(12)  we  obtain  an  equivalent  problem  to  (5), (6), (7)  of  the 
form: 


OJ)  Mu  l cTF*?  . [ (-1)cTF*2. 

♦,0  A 11 

subject  to 

(14)  0 1 1 dt  * dj 

I * n(0)  < j[|h  F. |z.<  g ‘ *l(0)  where  n * rank  (F) 

1*1 

The  system  (14),  (IS)  will  be  inconsistent  iff  either 
06)  t*  « n(0) 


or 

(17)  n(d*  • i~)  < l 

and  we  have 


Laima  1:  Necessary  and  sufficient  conditions  for  problem  (13) , (14) , (15) 

to  be  feasible  is  that  both  n(0)  i g * and  n(d*  - d~)  i g’. 

Proof.  Follows  directly  by  substituting  the  upper  and  lower  bounds  of  z in  the 
additional  constraint  (15). 

Let  ^ denote  the  ratio  of  the  coefficients  of  z A in  (13)  and  (15), 
then  Yj  can  be  written  as: 


» 


(18) 


i c {♦} 

i c (-) 

i e (0 ) and  cTF*  ^ 0 
i e (0 ) and  c*F*  < 0 


i 
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where  M is  a positive  nunber  of  dominating  magnitude. 


(At  this  jxmit  we  may  write  down  the 

optima1  values  of  z.  for 

1 L (0>  - > Z-  = 

' d+  - 
1 

d!  if  cTF#  > 0 

: l 

0 

k 

if  cTF#  < 0)  . 

l — 

Let 

r 

i e {♦,0> 

(19) 

H- 

T J 

i e {-} 

■cFi 

uiuwr  • ♦ , 0}  = { ♦ \ (j  { (1 } , 

and  let  (y  Ci)}  denote  the  set  of  Yj's  arranged  in  decreasing  order,  with 

(j)  * (k)  also,  if  y^  * Y^  but  the  original  in^ex  of  (j)  is  less  than  the  original 
index  of  (k)  .where  (i)  is  the  i*h  index  order  of  (r^h 

i ft 


Ji>  i 

p = 

Max  integer 

3 “(i)  - U; 

p = 

0 if 

rjtri)  < °* 

:i  i 

pi  = 

Least  integer 

Pi 

3 sP  = r 

1 i=l 

ihVa,i 

(d(«-  d'(i)’ > 

P2  = 

Ps 

l.east  integer  3 Sp.,  - 1“ 

“ i=l 

|hV(i)| 

(d(ij-  <w  > 

n 


|hTpri  \ I d'  ,)  f g+  - n(0),  set  P,  - n + 1. 


!f  Sn  - 7^  u(i)'  r & JV-U  *2 


Ii  ~ g - n(0),  set  Pj  * n ♦ 1.  If  Sn  < g - n(0),  the  problem  is  inconsistent. 

Remark  l 

Note  that  since  the  denominator  of  Yj  is  always  positive,  the  elements 
in  (y^)}  are  oidered  in  such  a way  that  those  with  positive  values  of  ct^ 
appear  first,  and  those  with  negative  values  later. 


Remark  2 


The  assunption  c*F^  f 0 can  always  be  made  if  desired  since,  if  not, 
a perturbation,  essentially  equivalent  to  that  introduced  by  Chames[2j  for 

the  simplex  method,  can  be  performed.  The  perturbed  problem  is  obtained 

T # 

from  the  original  problem  by  replacing  the  vector  c F by  the  perturbed  vector 
(23)  (c£). 


V * r1 
C r.  + c 


where  c is  a sufficiently  small  and  positive  number,  so  that  optimal  solutions 
to  the  perturbed  problem  are  optimal  solutions  of  the  original  problem 


and  (c#)  7*  0 (i  * l,...,n). 


The  optimal  solution  2*  to  (13) , (14) , (15)  can  now  be  written  for  the 
different  cases  in  the  following  manner: 

Case  I:  For  P < Pj  < the  optimal  solution  is: 


d - d,. 


(i)  (i) 


(24) 


•v* 


Pi-1 


1 < i < P,-l 

- - i 


‘w  ’ 1 Vi)  * Ig'  * "(0)  -jilhTF'(i)l(da)-d‘(i))]/|hTF,(Pi)|  "j . 


lTb# 


Pj+1  < i < n 


Case  II:  For  < P the  optimal  solution  is: 


(25) 


!(i)  " 1 


/*( i)  ‘ dCi) 


1 < i < P,-l 
- - L 


P?*l 


(p2)  ■ • n(0)  •iI1lhV(i)l(da)'  d(i))J/lhTF(P2)i  i - p. 


Tr# 


P?*l  s i < n 


l 


1 


? i 


Case  III:  For  < P < Pj  the  optimal  solution  is: 


1 < i < P 


P*1  < i < n 

• m 


Remark:  In  cases  I and  II  above,  the  denominator  of  the  expression  defining  0 

will  not  be  zero.  The  reason  for  this  is  as  follows: 


— -e  I:  h F(P1)  = 0 Px  = 1,  since  zero  is  a lower  bound  for 

p 

Spi  = ii1  F(i)^d(i)"  d(i)^  3,1(1  pi  is  the  least  integer  for  which  Spx  > g‘  - 

T # 

But  h F(I)  - 0 also  implies  that  c(1)  > 0,  from  the  definition  of  y i.e. 

This  contradicts  P < Pj  for  case  I. 

case  If:  h i-^p^  * 0 P2  ~ 1 (as  alx)Ve) 

Sp2  = 0. 

But  g - n(0)>  0 is  a necessary  condition  for  consistency  of  (14),  (15).  Hence 
SP2  ^ 8 * n(0).  This  contradicts  the  definition  of  P^. 


The  corresponding  optimal  values  of 


*(i)  for  the  three  cases  are: 


(i)e  {+,0} 

1 5 i < Pj-1 

(i)e  {-} 

1 < i < P,-l 

(i)c  {♦} 

-h 

(i)c  {-} 

i-Pi 

Ci)e  {*,0} 

^l+l  < i < n 

(*)«  {-> 

Pl*l s i i n 

n(0). 

P > l. 
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Case  II: 


(28) 


(i) 


Case  III: 


(i) 


% 

(i)«  (+,0) 

1 < 1 < P2-l 

d(i) 

(i)c  (-) 

1 < i < P2-l 

V2)’d(P2) 

(i)e  {♦} 

i-P2 

d(P2)'V2) 

(i)e  {-)  ■ 

i - P2 

d'(i) 

(i)e  {♦,0} 

P2*l  s i i n 

d(i) 

(i)c  {-) 

P2*l  < i < n 

d(i) 

(i)  e {+,0} 

1 < i < P 

d"> 

(i)  e {-} 

1 < i < P 

dU) 

(i)  e {+,0} 

P + 1 £ i in 

(i)  e (-} 


P + 1 < i < n 


Remark 


Observe  that  fl,n  d 


(pl)+  d(pi)  ^ iS  ^JldePenc^ent  °£  d(p  ) 


and  d 


(pi)  ‘ 


To  show  this  we  set 
(30)  . 9- 


CPX)  * “'(P/  '8'  - l *VA  ‘j1lhTF(i)l(d(i)'d’(i)) 


Pr1 

TeM  r i vTr^ 


* d(P1)hT’:(P1)1/  hTF(Pj) 
Since  e {+,0}  we  obtain  from  (30)  that 

(31)  V/V/'8':  J 01hMd'i1hTiK 

»<?!>  Prl 


-.IIhV(i)ICd*(i)-  a-(l)),A¥(Pi) 


which  is  independent  of  d'rv  •.  and  drD  % . 


",  l.> 


^ 4. 

In  the  same  way  we  get  that  d^p  ^ - ®q>  ^ is  independent  of  d^p^  and  d^j 
Old  that  9-^  d'py  and  d^}  - 9^,  are  Independent  of  d^j  and  d'^ . 
For  each  of  the  three  cases  above  the  optimal  solutions  x*  for  (5), 


(6), (7) are  given  by 
(32) 

where  y is  arbitrary. 
Example  1 
Solve 


X*  - Ff z*+(i  - F#F)v 


Hex  *j  ♦ 2*2 
subject  to 
-9  < -3x1  ♦ x2  < 9 
0 < x2  < 8 

2 s xx  ♦ x2  < 6 


In  this  case  we  choose  F 


• r.  ;)• 


We  have 


F*  * F'1  * ^ 

hT  - (1.1) 


-1/3  1/3 


cV1  - (-1/3,  7/3) 


hV1  - (-1/3,  4/3) 


Thus  the  problem  is  transformed  into 

Mgx  -l/3z^  ♦ 7/3z2 
subject  to 


■9  < z. 


?2  < » 


2 < -1/3ZJ  ♦ 4/3tj  < 6 


Now  let 
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'N*  - 

*1  " ^ " Z1 


l2  1 2 


The  problem  to  be  solved  is 


Max  -1/3(9  - Zy)  * 7/3z2  * Max  l/3zx  ♦ 7/3z2 
subject  to 

A. 

< 18 


_ % 
0 < Zj 

0 < 


-A#  A 

z2  < 8 


We  have 


Y1  * I7J  " 1 


. 5 < l/ZZy  ♦ 4/3z2  < S 


Y - 111  • 7/4 
Z 4/3 


{Y(i)}  " {7/4»1} 


Thus 


*i 


I |hTF(J)l  (d*}-  d'  ) - 4/3  • 8 - 10  | > 9 
i-1 

P ■ 2,  and 


-3/4 
27/4  / 


Theorem  1 

If  feasible,  an  optimal  solution  to  problem  (5) ,(6), (7)  is  given 
for  the  three  cases  I,  II,  III  by  (27), (28), (29)  respectively. 


A 
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I 

£ 

Proof:  Follows  directly  from  the  way  was  constructed, 

Remark  4: 

I 

If  (5), (6), (7)  is  feasible  it  is  always  possible  to  delete  a con- 

L 

straint  from  the  set  of  constraints  (6),  (7)  without  affecting  the  optimal 

p 

[ solution.  This  follows  from  Remark  3 ii  which  we  showed  that  in  Case  I, 

r 

z is  indepe  ndent  of  (Pj),  thus  the  (Pj)  constraint  of  F can  be  deleted. 

• 

j In  Case  II  constraint  F can  be  deleted,  and  in  Case  III ^ the  additional 

t ^ 

constraint  (7)  can  be  deleted.  In  each  of  the  above  cases  deletion  of  the 

* 

j suitable  constraint  will  not  change  the  optimal  solution. 

t 

1 


2 An  Algorithn  for  Solving  General  (IP)  Problems 

In  this  section  an  algorithm  for  solving  (IP)  problems  will  be 
described.  The  (IP)  problem  to  be  solved  is: 

(33)  Max  cTx 

subject  to 

(34)  b"  s Ax  < b* 

where  b’  < b+f  cJLN(A)  and  A e . If  m ■ r then  an  optimal  solution  to 
(33) , (34)  is  given  by  (29) . 

For  x c R1**  let 

(35)  I(x)  ■ (i:  (Ax)a  > b*  or  (Ax)A  < bT)  . 

Let  H(0)  be  any  set  of  r indices  from  (l,...,ra)  such  that  (a^:  i e H(0)) 
is  a linearly  independent  set  of  rows  of  A,  and  let 

(36)  i(0)  e (1 ra)  i(0)  t H(0) 

Proceed  to  the  first  iteration: 

Iteration  v > 1:  Apply  the  method  in  I to  determine  infeasibility  or 

to  find  an  optimal  solution  x*(v)  to  the  following  problem: 


r, 
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Problem  v: 


(37) 

Max  cTx 

subject  to 

(38) 

b[  < aAx  < b?  i e H(v-l) 

(39) 

bi(v-l)  - ai(v-l)x  - bi(v-l) 

If  problem  v is  infeasible  then  (IP)  is  infeasible;  stop. 

Alternately,  problem  (33)  (34)  can  be  "regularized"  a la  Chames-Cooper 
so  that  infeasibility  is  determined  at  optimal  solution  if  that  is  the  case. 

If  I(x*(v))  ■ t then  x*(v)  is  an  optimal  solution  of  (IP);  stop. 

Otherwise  let 

(40)  i(v)  c I(x*(v)) 

For  Case  I we  choose: 

(41)  H(v)  - H(v-l)  U (iCv-DJ/tPjCv)} 
where  / denotes  deletion. 

For  Case  II 

(42)  H(v)  - H(v-l)  U (i(v-l)  >/{P2(v)) 

(For  the  definition  of  P^(v)  and  P2(v)  see  (21), (22)).  For  Case  III, 
run  II, 'v)  = H(v  - 1).,  . 

j^ocoo.1  to  iteration  v + 1. 

Notice  that  from  the  definition  of  P^  and  P2  the  matrix  with  rows 
a^,  i c H(v)  is  of  full  row  rank.  Moreover,  since  the  coefficient  matrix 
in  iteration  v differs  from  the  coefficient  matrix  in  iteration  v-1  by  only 
one  row,  the  product  form  of  the  inverse  (e.g.  Chames  and  Cooper  J3] ) may 
be  used  to  compute  the  new  inverse. 

Theorem  2 

The  algorithm  described  above  terminates  in  a finite  mmber  of 
steps  either  with  the  conclusion  that  (IP)  is  infeasible,  or  with  an  optimal 
solution  to  (IP). 


V 'rt*'yTV4  -•  v 
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Proof:  Since  in  iteration  v the  algorithm  solves  a full  row  rank  problem 
with  one  additional  constraint  which  was  not  satisfied  by  the  optimal 
solution  x*(v-l),  it  follows  from  Remark  2 ttet 

(43)  cTx*(v-l)  > cTx*(v) 

which  assures  the  finiteness. 

The  optimality  is  assured  since  at  each  step  of  the  algorithm  we 

i 

find  an  optimal  solution  to  a less  restricted  problem. 


Example  2 

Solve 

Max 

X1  4 ^2 

subject  to 

0 5 

X1 

6 

• 

0 < 

8 

.2  < 

X1  4 *2  -< 

6 

-9  < - 

3xx  * *2  < 

9 

Iteration  0: 

• 

H(0)  - (1,2) 

i(0)  - J 

Iteration  1: 

Solve 

• 

Max  Xj  ♦ 2x2 
subject  to 
0 f Xj  s 6 
0 5 Xj  5 8 
2 s Xj  ♦ X2  s 6 

Px  - 1,  P2  - 1,  P - 2 

which  is  the  optimal  solution  since  I(x*(l))  * 0. 


(y(i))  - (2,1) 
From  (24)  x*(l)  «pj 
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Example  3 [ 5] 

Solve 


Iteration  0: 


Iteration  1: 
Solve 


(T(i)>- 

Thus 

Iteration  2; 

Solve 


Max  x,  ♦ 2x2 
subject  to 

0 < Xj  <6 

0 < x2  i 8 

•9  < -SXj  ♦ Xj  < 9 

2 < Xj  ♦ X2  < 6 


H(0)  - {1,2}  i(0)  - 3 


Max  -Zj  ♦ 2z2 

subject  to 

0 < _ < 6 

0 < \2  < 8 

*9  < &,  ♦ 2,  < 27 

1 Z - 


,-1/3)  P - 1 Pj  - 2 P2  » 2 

, and  I(x*(l))  - (4) 

Max  *2  + Zx2 

subject  to 
-9  < -3x,  ♦ x,  < 9 
0 < x2  < 8 

2 < X2  ♦ x2  < 6 
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Uuis  substituting  F(2)x  ■ Z we  obtain 

Max  -l/Stj  ♦ 7/3z2 
subject  to 

-9  < tx  <9 

0 < 

2 < -1/3*1  4V3z2  1 6 

we  obtain: 

Max  +\/l\  * 7/3z2 
subject  to 

& < 18 


substituting  \ 


z2  < 8 


0 < Zj 

0 < 

5 < l/l\  +4/3z2  < 9 


i2  < 8 


and  I(x£)  * CD* 
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Iteration  3: 


Solve 


Max  x i ♦ 2x2 


F(3)  - 


♦ 1/4z2 
1/Atx  * 3/4z2 - 


Mu  -l/4ix  ♦ 1/4z2  ♦ 2/4IJ  ♦ 6/4*2  - M“  l/4ll  * 7/4 z2 
subject  to 

z,  5 9 


*2  5 6 


substituting  1 


-9  s tx 
0 < 

0 < -l/4z1  ♦ 1/4z2  < 6 
we  obtain 

Max  1/4(9-^)  ♦ 7/4z2 


0 < 
0 < 


% 


< IS 
?2  < 6 


{T(i)}*  (7,-D 


9/4  < l/4zx  ♦ 1/4z2  <6*9/4 
P - 1 


Pi -2 


P2-2 
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I (x* (3) ) - 0-  ,nai»  x*(3)  is  optimal  for  (IP). 

Observe  that  the' lumber  of  iterations  needed  in  order  to  solve 
Example  2 was  strongly  influenced  by  the  choice  of  index  frqm  I(v). 
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